Steady state solutions of hydrodynamic traffic models 
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We investigate steady state solutions of hydrodynamic traffic models in the absence of any intrinsic 
inhomogeneity on roads such as on-ramps. ft is shown that typical hydrodynamic models possess 
seven different types of inhomogeneous steady state solutions. The seven solutions include those 
that have been reported previously only for microscopic models. The characteristic properties of 
wide jam such as moving velocity of its spatiotemporal pattern and /or out-flux from wide jam are 
shown to be uniquely determined and thus independent of initial conditions of dynamic evolution. 
Topological considerations suggest that all of the solutions should be common to a wide class of 
traffic models. The results are discussed in connection with the universality conjecture for traffic 
models. Also the prevalence of the limit-cycle solution in a recent study of a microscopic model is 
explained in this approach. 

PACS numbers: 05.45.-a, 89.40.-a, 45.70.Vn, 05.20.Dd 
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I. INTRODUCTION 

Vehicles on roads interact with other and various traf- 
fic phenomena can be regarded as collective behaviors of 
interacting vehicles Q, Q Q - Analyses of highway traf- 
fic data revealed that there exist qualitatively different 
states of traffic flow 0, IE IE • Transition between 
different traffic states was also reported Some em- 
pirical findings, such as the so-called synchronized traf- 
fic flow phase have become subjects of considerable 
concern and ignited intense theoretical investigations of 
traffic flow. Comprehensive reviews can be found for ex- 
ample in Ref. 

Many traffic ffow models have be en p ut forward and 
analyzed |li|ll|ll|ll|ll|ll|lllHlIlE|2E|2ll. 

Quite often, the analysis aims to find out possible steady 
states, or dynamic phases, of the models and investigate 
their properties. A criterion for good and reliable traffic 
models is good agreement between steady states of mod- 
els and traffic states revealed via real traffic data analysis. 
Thus an important first step in the analysis of traffic flow 
models is to find out all possible steady states. 

In this paper, we consider hydrodynamic models for 
traffic without bottlenecks such as ramps and present a 
systematic search for their steady state solutions which 
are time-independent in a proper moving reference frame, 

P{x,t) = P{X + Vgt), V{x,t) ^V{X + Vgt), (1) 

where p and v are density and velocity fields, respectively, 
and —Vg is the constant velocity of the moving reference 
frame with respect to the stationary reference frame. 
Two well-known steady solutions of this type are free 
flow and traffic jam solutions IHHIp. Surprisingly, 
we flnd that hydrodynamic models possess not only these 
two but also several other steady state solutions. Some 
of the newly recognized steady state solutions, including 
limit cycle solutions, have been reported previously only 
for microscopic traffic models 0,0, |2^ and not been re- 
ported for hydrodynamic models, which led to the wide 



spread view that free flow and traffic jam are the only 
possible steady state solutions of hydrodynamic models 
in the absence of bottlenecks such as ramps. Our result 
shows that such view is incorrect and that the physics 
contained in hydrodynamic models may be much richer 
than previously recognized. 

In Sec. ^ we flrst review the mapping to the single 
particle motion and introduce the concept of the flow 
diagram in the single particle phase space. It is demon- 
strated in Sec. IIIII that the flow diagram can have vari- 
ous topologically different structures, which are directly 
linked to the existence of certain types of steady state 
solutions CSec. IIV|I. Section E] discusses implications of 
our results. Section ?? concludes the paper. 



II. MAPPING TO SINGLE PARTICLE MOTION 

We consider a hydrodynamic model that consists of 
the following two equations, an equation for local vehicle 
number conservation, 
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where Vop {p~^) is the so-called optimal velocity function 
and the coefficients {TZ for the relaxation term, A for 
the anticipation term, and T) for the diffusion term) are 
positive deffnite and depend in general on the density and 
velocity fields, i.e., TZ{p,v), A{p,v), and 'D{p,v) > 0. 

To find out steady state solutions of the type in Eq. Q), 
it is useful to map the problem into a single particle mo- 
tion problem by using the method in Refs. 22, 24]. For 
the mapping, one flrst integrates out Eq. ^ . The result- 
ing constant of motion, 
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relates two dynamic fields p and w, and can be used to 
reduce the number of independent dynamic fields from 
two to one (we choose v in this work). Then Eq. (PJ 
can be transformed into an ordinary differential equation 
for the single dynamic field v that depends on a single 
parameter z = x -\- Vgt, 
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— +C{v;vg,qg)— +T{v;vg,qg) = 0, 
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Here the field p in the arguments of 7?., A, P, and Vop 
should be understood as qg/ {v + Vg). Then the search for 
steady state solutions reduces to the analysis of Eq. 
under the physically meaningful boundary condition that 
solutions should be bounded as z ^ ±oo. 

To gain insights into implications of Eq. jSJl, it is use- 
ful to make an analogy with a classical mechanics of a 
particle by regarding z as a time variable and w as a 
coordinate of a particle with unit mass moving in a one- 
dimensional system. Then Eq. ||SJ) describes the time 
evolution of a particle under the influence of a poten- 
tial energy, U{v;Vg,qg) = / dv' !F{v']Vg,qg), and of a 
damping force with the coordinate-dependent damping 
coefficient C{v; Vg,qg). For a physical choice of Vop{p~^), 
which decreases as p increases, goes to zero for large p, 
and saturates at a finite value for small p, the potential 
energy U becomes camelback- shaped for wide ranges of 
Vg and qg (solid curve in Fig. ^ . Thus the potential en- 
ergy profile is of very typical shape. Below we focus on 
U of this shape only and ignore the possibility of more 
exotic shaped U's, such as U's with three peaks, since we 
do not know of any reason to expect such exotic possi- 
bilities. Furthermore we will not address the trivial case 
occurring when the range of Vg and qg allows less than 
two peaks in U. 

An unusual feature in this mechanical analogy is that 
the damping coefficient C is not necessarily positive and 
in those ranges of v, where C is negative, the particle may 
gain energy due to the damping. The possibility of the 
negative damping is crucial for the existence of certain 
steady state solutions presented in the next section. 

Before we close this section, we remark that the stabil- 
ity of a solution in the single particle problem should 
not be identified with the stability in the real traffic 
problem. We demonstrate this point with trivial z- 
independent solutions. For the camel-back-shaped po- 
tential in Fig.^ extremal points of the potential become 
solutions. Thus there are three z-independent solutions, 
V = Vi (i = 0, 1, 2 and < vi < vq < where vq is the 
coordinate of the local minimum and vi, V2 are coordi- 
nates of the two local maxima or saddle points. All w^'s 
satisfy Vop{{v + Vg)/qg) = v and thus depend on Vg and 
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FIG. 1: A typical shape of the potential energy U{v;Vg,qg) 
(solid line). Three examples of the particle motion are also 
shown (dotted, dashed, dash-dotted lines). 



qg. These z-independent solutions correspond to the ho- 
mogeneous traffic states v{x,t) — wo,i,2- From the shape 
of the potential energy, it is clear that in the single par- 
ticle problem, the two solutions v{z) = vi^2 are unstable 
with respect to small deviations and the other solution 
v{z) = vq is stable (if C is positive near vq). For the 
real traffic problem [described by Eqs. ^ and Q], how- 
ever, the solutions v{x,t) = Vi^2 are (usually) stable with 
respect to small deviations, and v{x,t) — vq is linearly 
unstable. 



III. FLOW DIAGRAMS 

Besides the trivial static (z-independent) solutions, 
there also exist dynamic (z-dependent) solutions, some 
examples of which are shown in Fig. ^ (dotted, dashed, 
and dash-dotted lines). In the language of traffic flow, 
these dynamic solutions correspond to steady but inho- 
mogeneous traffic flow states. For a systematic study 
of dynamic solutions, it is useful to introduce a two- 
dimensional phase space {v,w = dv/dz), where each tra- 
jectory in the phase space corresponds to a dynamic so- 
lution. Then finding all solutions for given Vg and qg is 
equivalent to constructing a flow diagram in the phase 
space for the given Vg and qg. 

To gain an insight into flow diagram structures, it is 
useful to examine the flow near the three fixed points 
{v = Vi,w = 0)'s. Figure 121 show flows near the three 
flxed points. The fixed points (i'i.2,0) are saddle points 
regardless of the sign of C near v — vi^2 since the damping 
force alone [the second term in Eq. Q] can not reverse 
the direction of the particle motion (that is, the sign of 
dv/dz) even if the damping coefficient is negative. On the 
other hand, the fixed point {vo,0) is a stable (unstable) 
fixed point of the flow diagram when the C is positive 
(negative) near v = vq {C is assumed to be positive in 




FIG. 2: A schematic diagram of the fiow near the three fixed 
points (110,1,2,0). C is assumed to be positive near v — vq in 
this figure. 



Fig. I2J). Flows running out of or into the fixed points 
can be mutually interconnected and the way they are 
interconnected will in general depend on the values of Vg 
and Qg and affect the structure of the flow diagram. 

To make our discussion concrete, we choose here a par- 
ticular hydrodynamic model. We choose the coefficients 
in Eq. (O as follows: 



7^ = A, A 
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where 77 = dVop/d{p^^). As demonstrated in Ref. [2^, 
this choice provides a macro-micro link between the hy- 
drodynamic model [Eqs.^J and (PJ] and the microscopic 
optimal velocity model Q, 



y„ = A [14p(A?/„) - y„ 



(8) 



where ?/„ (t) represents the coordinate of the n-th vehicle 
in a one-dimensional road and Ay„ is the distance to the 
preceding vehicle y„+i — j/„. But we remark that as far 
as steady state solutions are concerned, the choice {T)) 
is just one particular option and most results presented 
below are not sensitively dependent on it. Results which 
are dependent on the choice will be stated so. For the 
parameters, values in Ref. |23| are used: A = 2 (sec~^). 
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33.6 (m/sec), y^d = 25.0 (m), ywidth = 23.3 (m), 
and Cicf ~ 0.913. For this Crcf, the maximum value of 
Vop is [(H-0.913)/2]i;inag, which is slightly different from 

For this hydrodynamic model, the resulting flow dia- 

and Qg in Fig. O 
and (7g, trajectories depart- 
ing from the two saddle points (wi,2, 0) behave in different 
ways and thus the flow diagrams acquire topologically dif- 
ferent structures. Since the structure of the flow diagram 
is closely linked to characters of nonhomogeneous steady 



gram is shown for various values of Vg ana Qg m 
Note that depending on Vg and qg 




FIG. 3: Schematic drawings of the fiow diagram in the phase 
space {v,w). Depending on the values of Vg and qg, the flows 
running out of or into the fixed points are connected in differ- 
ent ways and as a result, the fiow diagram acquires different 
topological structures. In each flow diagram, the fixed points 
(wo, 1,2,0) are marked by black circles. Emphasis is given to 
the fiows running out of or into the fixed points in the range 

Vl < V < V2. 



state solutions, it will be meaningful to divide the plane 
(vg,qg) according to the flow diagram structures, which 
is given in Fig. ^ The Vg-qg plane is divided into six re- 
gions (region I, II, • • ■, VI). The flow diagram within each 
region is labeled accordingly in Fig. O On the boundary 
between two neighboring regions, (eg. boundary Biji 
between the region I and II) the flow diagram acquires 
structures topologically different from those within the 
regions, and on the special point (t;*,g*), where all six 
boundaries join together, the flow diagram acquires a spe- 
cial structure still different from all others. Steady state 
solutions contained in the flow diagrams are presented in 
the next section. 



IV. STEADY STATE SOLUTIONS 

Out of all flow trajectories contained in the flow di- 
agrams (Fig. |3Jl, only those trajectories that remain 
bounded both for z — > 00 and —00 constitute physically 
meaningful steady state solutions. Below we focus on 
those bounded trajectories. 



A. Saddle— Minimum solution 

For each {vg, qq) in the region I and II, there exist a sin- 
gle trajectory starting from the saddle point (w2,0) and 
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FIG. 4: The division of the parameter space {vg,qg) based on 
the topological structure of the flow diagram. 
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FIG. 5; Inhomogeneous steady state solutions of a hydrody- 
namic model. 

converging to the potential minimum point (wo,0) [see 
Figs. I2l,ll)]. This trajectory represents the steady state 
solution in Fig. [Sfa). Similarly for each (vg,qg) in the 
region II and III, there exist a single trajectory starting 
from the other saddle point {vi , 0) and converging to the 
potential minimum point {vo,0) [see Figs. |3JlI, III)]. The 
steady state solution for this trajectory is similar to that 
in Fig. |3Ja) except that the spatial profile approaches vi 
instead of 1)2 as z ^ — cx). We call this type of steady 
state solutions saddle-minimum solutions. We remark 
that the oscillation near vq may or may not appear de- 
pending on the parameter choice, which can be easily un- 
derstood in the particle analogy; If the particle motion 
near vo is underdamped/overdamped, the convergence 
to Vo in the saddle- minimum solution is oscillatory/non- 
oscillatory. Expansion of Eq. (5) near {vo,0) shows that 
the motion near vq is overdamped if A > 2 and under- 
damped if A < 2 

A study 25] on the optimal velocity model revealed 
that the microscopic model possesses so-called oscilla- 



tory traveling wave solutions and monotonic wave solu- 
tion. These solutions are identical in their character to 
the saddle-minimum solution with the underdamped and 
overdamped convergence to vo ■ 

B. Limit-cycle solution 

For each {vg,qg) in the region IV, V, or VI, there ex- 
ists a single trajectory, which encircles the potential min- 
imum point (wo,0) and makes a loop [see dashed curves 
in Figs.|2IIV,V,VI)]. In the language of nonlinear dynam- 
ics, this type of flow is usually referred to as a limit-cycle. 
For traffic flow language, the limit-cycle corresponds to 
a steady state solution with periodic wave [Fig. IHb)]. 
The "wavelength" of the limit-cycle solution increases 
and approaches infinite as {vg,qg) approaches {v*g,q*g). 
We remark that for the existence of the limit-cycle so- 
lution, it is crucial to have sign alternation of C with v. 
If C were always positive (negative), the single particle 
motion analogy indicates that the total energy of the par- 
ticle would monotonically decrease (increase) with z, and 
thus the trajectory would be attracted to (repelled away 
from) destroying the limit-cycle. 

Note that since there is a limit-cycle solution for each 
(wg, qg) within the region IV, V, and VI, there is an infi- 
nite number of limit-cycle solutions, each with different 
{vg, qq). This feature is very similar to the report of many 
stable nonhornoeeneous states in a revised car-following 
model in Ref. [Tg ■ To our knowledge, it has not been re- 
alized previously that hydrodynamic models also possess 
infinitely many limit-cycle solutions. 

C. Limit-cycle— Minimum solution 

The limit-cycle solutions are inevitably accompanied 
by still different types of solutions. According to the 
flow diagrams in Figs. |3lV,V,VI), all trajectories inside 
the limit-cycle approach the minimum point (wo,0) as 
z — > 00, and the limit-cycle as z — ^ —00. We call this 
type of steady state solutions limit-cycle-minimum solu- 
tion. Its profile is shown in Fig. [SJc). Since there are 
infinitely many limit-cycle solutions and the limit-cycle- 
minimum solution is possible whenever the limit-cycle 
is possible, there also exist infinitely many limit-cycle- 
minimum solutions. 

This solution is closely related with an interesting phe- 
nomenon reported in previous studies of hydrodynamic 
models with on-ramps 14,„ J_^] : When a linearly unstable 
but convectively stable homogeneous fiow is generated 
in the upstream part of an on-ramp, oscillatory flow is 
spontaneously generated from the homogeneous region 
and propagate in the upstream direction. The relation 
between this observation and the limit-cycle-minimum 
solution can be established by noting that the convec- 
tively stable homogeneous flow region corresponds to the 
trivial solution v{z) = vq and the spontaneously gener- 
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ated oscillatory flow to the limit-cycle solution. It is then 
clear that the phenomenon in Refs. 0,0] is nothing but 
a manifestation of the limit-cycle-minimum solution. Re- 
cent studies of a microscopic car-following model 16] also 
reported oscillatory flow growing out of a linearly unsta- 
ble homogeneous region. 

D. Limit-cycle— Saddle solution 

Each limit-cycle solution accompanies still another 
type of solutions. For each {vg ,(?£,) in the region IV and V, 
there exists a single trajectory converging to the saddle 
point {vi ,0) as 2; — > oo and approaching the limit-cycle 
as z ^ —00 [Figs.|3|[IV,V)]. Similarly for each {vg,qg) in 
the region V and VI, there exist a single trajectory con- 
verging to the other saddle point {v2 ,0) as z — > 00 and 
approaching the limit-cycle as z —00 [Figs. EIVjVI)]. 
We call this type of solutions limit-cycle-saddle solution. 
Its profile in the z-space is shown in Fig. Hfd). Simi- 
larly to the limit-cycle-minimum solutions, there exist 
infinitely many limit-cycle-saddle solutions. 



E. Hetero saddle-saddle solution 

Still new types of bounded trajectories appear on the 
boundary B^j between the region i and j. On the bound- 
aries Biji and Biv,v, there is a trajectory running from 
the saddle point (wi, 0) and converging to the other sad- 
dle point (v2,0) [Figs.|2lBiji),(Biv,v)]- This trajectory 
amounts to the kink solution in Fig. [SJe). Also on the 
boundaries Bujh and Bv,vi, there is a trajectory run- 
ning from the second saddle point {v2 , 0) and converging 
to the first saddle point (wi,0) [Figs. |3|[Bnjii),(Bv,vi)]- 
Since the roles of {vi,0) and (w2,0) have been swapped, 
this trajectory amounts to the anti-kink solution [not 
shown in Fig. 0. In the nonlinear dynamics language, 
these kinds of trajectories connecting two different fixed 
points are called heteroclinic orbits. In this paper, we 
will name this type of solutions hetero saddle-saddle so- 
lutions. To be more specific, we will call the first (sec- 
ond) type of the hetero saddle-saddle solutions "upper" 
("lower") hetero saddle-saddle solutions since the tra- 
jectories appear in the upper (lower) half of the flow 
diagram. Note that there exist infinitely many hetero 
saddle-saddle solutions since this solution is allowed for 
each point (vg,qg) on the abovementioned boundaries. 

F. Homo saddle— saddle solution 

For each {vg,qg) on Bvi,i, there exists a trajectory 
starting from the saddle point (u2, 0) and returning back 
to the same saddle point [Fig. IH^Bvij)]. This trajectory 
represents the narrow cluster solution [Fig. EJf)]. For 
each (vg,qg) on Bmjv, on the other hand, there exists 
a trajectory starting from the saddle point (wi,0) and 



returning back to {vi,0) [Fig.|3fBiiijv)]- This trajectory 
represents the narrow anti-cluster solution (not shown 
in Fig. |SJ). In the nonlinear dynamics language, these 
kinds of trajectories are called homoclinic orbits. In this 
paper, we will name this type of solutions homo saddle- 
saddle solutions. To be more specific, we will call the 
first (second) type of the homo saddle-saddle solutions 
"right" ("left") homo saddle-saddle solutions since the 
trajectories involve the saddle point on the right (left) 
half of the fiow diagram. 



G. Wide cluster solution 

Figure 01 shows that all six regions and six boundaries 
meet together at a single point {v* ,qg)- The flow diagram 
at this point has a special structure [see Fig.^v*,q*)] 
that allow smooth connection between flow diagrams 
of different topological structure in different regions or 
boundaries. 

This point is also special in the sense that the flow 
running out of the fixed point (^1,0) reaches the other 
fixed point (w2, 0) and returns back to (wi, 0). The steady 
state solution corresponding to this flow is shown in 
Fig. El^g). Note that since the particle dynamics be- 
comes infinitely slower as a trajectory approaches the 
fixed points (ui.2,0), the cluster size of the steady state 
solution is infinitely large. For this reason, we call this 
solution wide cluster solution. 

The wide cluster solution is a limiting case of various 
solutions mentioned above; If {v*, q*) is regarded as a lim- 
iting point of, for example, the boundary Bnijv or Bvi,i, 
the wide cluster solution is a homo saddle-saddle solu- 
tion with an infinite cluster size. If (v*, q*) is regarded as 
a limit point of the region VI, V, or VI, the wider cluster 
solution is a limit-cycle solution with an infinite period. 
Also if {v*,q*) is regarded as the point where the border- 
lines Bn.iii and Bi n join, the wide cluster solution is a 
combined object of the upper hetero saddle-saddle solu- 
tion (kink) and the lower hetero saddle-saddle solution 
(anti-kink). 

The fact that the wide cluster solution is possible only 
at the single point («*,(/*) implies that the wide cluster 
solution has the so-called universal characteristics; when 
a given initial state of traffic evolves into the wide cluster 
solution following the real traffic dynamics [Eqs. (|2Jl and 
(O], characteristics such outflow and the speed of the 
final traffic state are independent of the initial state. 

Analyses 0, of empirical traffic data revealed that 
various characteristics of the wide jam are indeed uni- 
versal. This empirical observation imposes a constraint 
on traffic models and the universality of the character- 
istics has been tested for traffic models. Although the 
universality has been verified for a number of traffic mod- 
els ^3 113 J the verification unfortunately has re- 
lied largely on repeated numerical simulations and thus 
the verification itself is also "empirical" in some sense. 
One exceptional case is the analysis in Ref. where 
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the singular perturbation theory is used. This analysis 
is however considerably model-dependent and thus it is 
not easy to perform the same analysis for other classes 
of traffic models. Contrarily, in our approach, the new 
discovery that the wide cluster solution can exist only at 
the single point {v*,q*) is quite meaningful since it ex- 
plains clearly and unambiguously why the characteristics 
of the wide cluster solution should be universal. 

We next address an interesting size dependence of the 
"wide cluster solution" reported in Ref. When nu- 
merical simulations are performed with periodic bound- 
ary conditions, the so-called universal characteristics are 
found [23, H^l to be not strictly universal but depend on 
the system size. This size dependence is not compatible 
with the statement of the universality given above. To 
resolve this conflict, we first note that periodic bound- 
ary conditions allow only those solutions whose periods 
are compatible with the imposed period. Then the wide 
cluster solution, whose period is infinite, can not be real- 
ized in such numerical simulations with finite system size, 
and the solutions, which are interpreted as the wide clus- 
ter solution in Ref. 22, 23], are, strictly speaking, limit- 
cycle solutions whose periods are compatible with the 
imposed periodic boundary conditions. Then by recall- 
ing that there are infinitely many limit-cycle solutions, 
the size dependence is a very natural consequence and 
the conflict is resolved. We remark that the dependence 
however should become weaker as the period becomes 
longer since the limit cycle solutions approach the wide 
cluster solution as their period becomes longer. 



V. DISCUSSION 
A. Roles of topology: Generality 

As demonstrated in previous sections, each steady 
state solution exists only in restricted regions of the Vg-qg 
space where the flow diagram acquires certain topolog- 
ical structures. This relation between the steady state 
solutions and the topological structure of the flow dia- 
gram is not a mere coincidence; all steady state solutions 
presented in the preceding section are guaranteed to ex- 
ist by the topologies of the flow diagrams. For example, 
when flows near the fixed point {vo,0) is attracted to- 
wards (wojO) while flows running out of the other fixed 
points (wi.2,0) are repelled away from (wo,0), as in the 
regions (IV,V,VI), there should exist the hmit-cycle solu- 
tion in those regions (Poincare-Bendixson theorem ^2l|)- 

This relation with the topology bears an interesting im- 
plication. As demonstrated in many branches of physics, 
physical objects, whose existence is closely related with 
a certain topological structure of systems, are not fragile 
and their existence does not depend on quantitative de- 
tails of the systems. Vortices in type II superconductors 
are a well-known example. It is then expected that the 
steady state solutions presented in the preceding section 
are not specific to the particular model examined but 



common to many versions of hydrodynamic models. For 
example, our results are not sensitive to the values of the 
parameters in Eq. 0. 

The generality due to the topology can be argued in 
the following way as well. Let us consider an infinitesimal 
change in the traffic model [in Eqs. Q and for exam- 
ple]. At (vg, Qg) located sufficiently interior of a region in 
Fig. 01 the topological structure of the flow diagram will 
not be affected and thus all steady state solutions, which 
are originally allowed in the region, are still allowed for 
the modified traffic model. On the other hand, at (vg, qg) 
located sufficiently close to a boundary in Fig. 01 the 
topological structure of the flow diagram may be modi- 
fied to a new structure. However the only possible new 
structure is the one at just across the boundary. The 
net effect then amounts to a mere shift of the boundary. 
Thus as far as the topology is concerned, effects of the 
infinitesimal change in the traffic model is no more than 
shifting the boundaries of the six regions by infinitesi- 
mal amounts, and the existence of steady state solutions 
in each region is not affected. Here we remark that de- 
spite the shifts, all six boundaries should still meet at a 
single point, whose coordinate may be slightly different 
from the original though. Otherwise, continu- 

ous conversion of one topological structure to another 
near (Vg^q*) is not possible. 

Lastly we discuss briefly one special kind of modifica- 
tions that shows negative-valued coefficient C [Eq. lO] 
near the fixed point {vq , 0) . In the hydrodynamic model 
proposed by Kerner et al. 0]j C is indeed negative when 
the parameters in the model are set to the values sug- 
gested in Ref. For this case, the flow is repelled 
away from (vq, 0). Then the limit-cycle solution does not 
appear any more in the regions IV, V, VI but appears 
instead in the regions I, II, HI. This shift is natural in 
view of the Pointcare-Bendixson theorem. However this 
negative C does not affect the very existence of the six 
regions because these are determined according to the 
trajectories departing from the two saddle point (wi,2,0) 
as mentioned already. Thus our approach focused on the 
bounded trajectories between the four limiting behaviors 
(one limit-cycle and three flxed points) is still valid. A 
more exotic kind of modifications are those, due to which 
C changes its sign multiple times with v so that multiple 
limit-cycle solutions exist for given Vg and qg in certain 
regions of parameter space {vg,qg). However this possi- 
bility seems to be very unlikely since it requires consider- 
able fluctuations of A [see Eq. ©], which is unphysical. 



B. Implications on universality conjecture 

A conjecture has been put forward by Herrmann and 
Kerner [29|; many traffic models with different mathe- 
matical structures may belong to the same "universal- 
ity" class in the sense that they predict same traffic phe- 
nomena. Although it is not clear yet to what extent the 
universality conjecture is valid, there are indications that 
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there indeed exists close relationship between some trafSc 
models. For example, macroscopic hydrodynamic mod- 
els have been derived from a microscopic car-following 
model via certain approximation methods '2^, l30l | , and 
good agreement between two types of traffic models has 
been demonstrated via numerical simulations |3ll |. In 
particular, an exact correspondence between two differ- 
ent types of microscopic models has been established [s^ . 
We note in passing that the hydrodynamic model derived 
in Ref. [30| however has an instability that arises from 
short-wave length fluctuations that is not present in the 
original car-following model. Such instability disappears 
when the short-wave length fluctuations are properly reg- 
ularized as in Ref. 26]. 

On the other hand, there exist reports which could 
not have been reconciled with the universality conjec- 
ture. For example, a recent study on a certain spe- 
cial type of car-following model reported the existence of 
many limit-cycle solutions (even without intrinsic inho- 
mogeneities on roads such as on-ramps). To our knowl- 
edge, the limit-cycle solution had not been reported for 
any other traffic models and thus it had been inferred 
that the limit-cycle solution might be specific to the spe- 
cial model, in clear contrast with the universality con- 
jecture. Our results however indicate that this inference 
is wrong and reopen a possibility that the special car- 
following model may also be closely related to other traf- 
fic models. An evidence for this will follow in the next 
subsection. 

By investigating topological structures of the flow di- 
agrams, we found from a single traffic model seven in- 
homogeneous steady state solutions. Although many of 
these solutions were already reported by earlier studies, 
earlier reports were scattered over various different traf- 
fic models and thus it was not clear whether a certain 
solution is specific to certain traffic models or generic 
to a wide class of models. Our analysis in the preced- 
ing subsection indicates that the seven inhomogeneous 
steady state solutions are generic to a wide class of traf- 
fic models. Our results are thus valuable in view of the 
universality conjecture. 

C. Prevalence of limit-cycle solution in Ref. |18|I 

As remarked above, most studies of trafhc models 
failed to capture the limit-cycle solution. On the other 

handj in the special car-following model studied in 

Ref. Iq , a wide class of initial conditions evolve to the 
limit-cycle solution. A question arises naturally: what 
is special about the model in Ref. 18]? The model is 
defined as follows. 





FIG. 6: The bold solid line is the effective optimal veloc- 
ity Vop due to the dynamics of Eq. (??). The dotted line 
on the other hand describes the constant of motion line 
in Eq. (|1J. The hatched region, where no dynamics is al- 
lowed, is established by the singular effect of the second 
term in Eq. (??). Three crossing points of the two lines de- 
termine the z-independent solutions of Eq. ©. Note that 



vehicles, Wpor is the permitted velocity. A: is a constant, 
and Ay^ = ijnT + p"^. Here T is the safety time gap. 
When interpreted in terms of the hydrodynamic model 
in Eq. JSJ, the first and third terms together define the 
effective optimal velocity V^p{p~^), which is shown in 
Fig. ??. The role of the second term is to strictly pre- 
vent the distance from being smaller than the minimum 
distance by establishing additional strong decelera- 
tion when a vehicle is faster than its preceding one and 
their separation approaches p^. 

In view of the steady state analysis given in previ- 
ous sections, the optimal velocity profile in Fig. ?? is 
very special: out of three solutions vi,V2,V3 of V^p{{v + 
Vg)/qg) — V, which amount to the extremal points of the 
potential , vi is dynamically forbidden since the ve- 
hicle spacing related to vi via Eq. is shorter than 
which is strongly prohibited by the second term in 
Eq. (??). Then steady state solutions, such as wide clus- 
ter solution, that assume the accessibility to {vi,0) are 
not allowed any more and the number of possible solu- 
tions reduces to 5 [saddle-minimum, limit-cycle, limit- 
cycle-minimum, limit-cycle-saddle, and homo saddle- 
saddle solutions]. Further reduction occurs when the 
periodic boundary condition is imposed as in Ref. |18| ]. 
Then homo saddle-saddle and limit-cycle solutions are 
the only possible solutions. Among these two, the for- 
mer is possible only when the average density satisfies 
P < (Pm^ + Tvpc,:)~^ because the unique limiting behav- 
ior of that solution, of which density converges to p2, will 
govern the average density and P2 < Pc- Thus for the den- 
sity range of {p^ + Ttipor)"^ < P < Pm, the limit-cycle 
solution is the only possible solution. Therefore the char- 
acteristics of the model in Ref. |I2j is explained within 
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the framework of the hydrodynamic approach. 

D. Stability of solutions 

The rigorous stabiHty analysis, whether each solution 
above is realized through the dynamics in Eqs. Q and 
0, goes beyond the scope of this work. In this section, 
we instead summarize what has been known and also 
discuss implications of existing results. The stability of 
the solution in Sec. II VI G is well-established ^,22\. For 
the solution in Sec. IIVI C, its relation with the oscilla- 
tory flow generated spontaneously out of a convectively 
stable homogeneous flow, which was reported in previ- 
ous numerical simulations of hydrodynamic models with 
on-ramps 0, 0|, seems to indicate that this solution 
can be stable. Also the solution in Sec. IIVI B have been 
maintained stably in our own numerical simulation of the 
hydrodynamic model that is derived from the following 
microscopic model via the mapping rule in Ref. |26l |: 

yn = \ [V,f (Ay„) - y„] + k _^ 9 {~Ay^) , 

(11) 

where V°p is the same one depicted in Fig. ??, fc is a con- 
stant, and 8(x) is the Heavyside function (1 for x > 
and for a; < 0) . Note that the idea of the prevalence of 
the limit-cycle solution discussed in Sec. C is simply 
reflected in Eq. (??). These observations suggest that at 



least some steady state solutions addressed in this work 
can be maintained stably. However more systematic anal- 
ysis is necessary to clarify the issue of the stability. 



VI. CONCLUSION 

Hydrodynamic traffic models are investigated by map- 
ping them to the problem of single particle motion. It 
is found that typical hydrodynamic models possess seven 
different types of inhomogeneous steady state solutions. 
Although these solutions were already reported by ear- 
lier studies, earlier reports were scattered over various 
different traffic models and it was not clear whether a 
certain solution is specific to certain traffic models only 
or generic to a wide class of models. Our result combined 
with the topology argument indicates that the seven in- 
homogeneous steady state solutions should be common 
to a wide class of traffic models. Also the origin of the 
universal characteristics for the wide cluster solution is 
clearly identified and the reason for the prevalence of the 
limit-cycle solution in a previous report [l^ is provided. 
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